Raman Spectra of Triplet Superconductor in Sr2Ru04 
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We study the Raman spectra of spin-triplet superconductors in Sr2Ru04. The p-wave and f-wave 
symmetries are considered. We show that there is the clapping mode with frequency of \/2A(T) 
and 1.02A(T) for p-wave and f-wave superconductors, respectively. This mode is visible as a huge 
resonance in the Big and B2g modes of Raman spectra. We discuss the details of the Raman spectra 
in these superconducting states. 
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I. INTRODUCTION 

The superconductivity in Sr2Ru04 was discovered in 
1994. El Shortly after the discovery of superconductiv- 
ity, a possible triplet p-wave superconductivity with the 
following order parameter was postulated by Rice and 
Sigrist.fl I 



A(k) = Ad(k x ± iky), 



(1) 



where A is the magnitude of the superconducting order 
parameter. Here d-vector is called spin vector perpendic- 
ular to the direction of the the spin associated with the 
condensed pair.JMJ Notice that this state is analogous to 
the A phase of ^le and there is a full gap on the Fermi 
surface. Subsequently, the triplet superconducting na- 
ture has been confirmed by the constancy of 17 Knight 
shift (spin susceptibility) through T c for the magnetic 
field parallel to the a — b plane [^| . The broken time rever- 
sal symmetry state have been also confirmed by the spon- 
taneous magnetic moment found in /iSR experiment ||. 
In general, the triplet superconductors have a variety of 
collective modes. The spin waves and the clapping mode 
with the order parameter of Eq. [j] were studied in Ref. 
J?], ||, [|. Further, the effect of the clapping mode on the 
sound wave was studied. || However, this coupling is very 
small to detect the existence of the clapping mode. 

In the meantime, the sample quality of Sr2Ru04 has 
been improved. The cleanest sample shows the transi- 
tion temperature close to the optimal T c — \.bK deduced 
from the T c dependence of the residual resistivity. [jl0| All 
these high quality samples exhibit the character of nodal 
superconductors; the T 2 behavior of the specific heat [p"l| , 
the T-linear dependence of superfluid density |l2|], the 
T 3 behavior of l/Tx in NMR||, the T 2 behavior 
of the ultrasonic attenuation Q, and the yflf depen- 
dence of the specific heat in a magnetic field [jll| at 
low temperature. These low temperature behaviors of 
the specific heat and superfluid density are consistent 
with a f-wave superconductor with nodes in the order 
parameter. |0| |j| [l7|, [jj, |l| 

On the other hand, the quasi-two dimensional system 
with strong paramagnon (the ferromagnetic spin fluc- 
tuation) favors the p-wave superconductor as given in 



Eq. [l].[^0| |T| In order to understand this situation, the 
nodal structure of the order parameter becomes of cru- 
cial importance. In a series of papers, Won and Maki 
have shown that the magneto-thermal conductivity for 
T <C v\[eH <C A(0) will provide the direct access to 
this question^, f23[ , where v — ^Jv a v c and v a . c are the 
anisotropic Fermi velocities. Such experiments have been 
carried out by Tanatar et a/[^4) and Izawa et al^Z^ . The 
result is consistent with the f-wave superconductor with 
horizontal nodes of Eq. ||, even though the thermal con- 
ductivity data cannot exclude a small admixture, a few 
% of p-wave order parameter, Eq. |l|. 



A(k) = Ad[(k x ± ik y ) cos (ck z )], 



(2) 



where c is the lattice constant along the c-axis. 

Parallel to this development, Zhitomirsky and Ricep(| 
have proposed an alternative model, multi-gap model for 
Sr2Ru04. As it is known that the Fermi surfaces in 
Sr2Ru04 consist of three different bands labeled by a-, 
/?-, and 7-bands||, it is also believed that the supercon- 
ductivity arises mainly in the 7-band. It was proposed 
that a full gap with p-wave, Eq. (1) exists in the active 
band 7, while line nodes with f-wave order parameter, Eq. 
^ develops in the a and j3 due to proximity effect. |26| 



A zr (k) = Ad[(k x ± ik y ) cos (ck z /2)}. 



(3) 



While this model could reproduce the specific heat data 
by Nishizaki et al [p"i"|| , and the magnetic penetration 
depth data by Bonalde et al|l^, |27|, they have not at- 
tempted to calculate the magneto-thermal conductivity 
which should be more revealing. Note that the f-wave su- 
perconductor associated with a and /3-bands in Ref. |2(J 
is similar to the Eq. ^|, but not the same. The model with 
A zr produced much larger (~ 30 times) cos(2<j)) term in 
the angular dependence of magneto-thermal conductiv- 
ity. The better test of the nodal position can be done 
through the angle (azimuthal angle of k from the a-axis) 
dependence of the thermal conductivity Jl7|, p8|. 

More recently, Deguchi et al observed a double transi- 
tion in the specific heat with the magnetic field near H C 2, 
and this experimental result was interpreted in terms of 
multigap model. |^| On the other hand, the behavior of 
the specific heat and the magnetic penetration depth for 
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FIG. 1: Order parameters with the cylindrical Fermi surface 
for (a) p-wave (b) f-wave of the Eq. g (c) f-wave of the Eq. |^ 



low temperature (T <C T c ) and the low field {H <C H C 2) 
appears to be consistent with the f-wave order parameter 
given in Eq. ^. At the moment, it is not clear whether the 
order parameter of Eq. ^, or the multigap model with the 
order parameters of Eq. [I] and Eq. yis adequate. This 
is the fundamental issue for Sr 2 Ru04. 

In this paper, we study the Raman spectra of spin 
triplet superconductors in Sr 2 Ru04. We consider three 
different order parameters given in Eq. |l|, 0, and || 
sketched in Fig. 1; (a) p-wave with a full gap, (b) f- 
wave with horizontal nodes, and (c) f-wave with node at 
the zone boundary. We identify the clapping mode[^ in 
these superconductors and consider its effect on Raman 
spectra. We show that (1) the Raman spectroscopy can 
detect the clapping mode, and (2) it can discriminate the 
multigap model from the single gap model. For the p- 
wave superconductor, the frequency of the clapping mode 
was found at \/2A(T)Q, and this mode exists as a sharp 
resonance in B\ g and Bi g modes of Raman spectra. For 
the f-wave superconductors (Eqs. || and ^|), the frequency 
of the clapping mode gets smaller than that of the p-wave, 
1.02 A as expected due to the existence of node, and this 
mode is also detectable in B\ g and B 2g modes of Raman 
spectra. While the weighted sum of the clapping mode 
contributions from the p-wave and the f-wave appears in 
the multigap model, the single gap model of Eq. || has 
the contribution from solely the clapping mode in the 
f-wave superconductor. 

The paper is organized as follows. The formalism to 
identify the clapping mode and its effect on the Raman 
spectra is summarized in the section II. The results of the 
Raman spectra of B\ g and B2 g modes for p-wave and f- 
wave superconductors are presented in the section III and 
IV, respectively. The conclusion and discussion will be 
followed in the section V. The details of computing the 
correlation functions are presented in Appendix. 



II. FORMALISM OF RAMAN SPECTRA AND 
CLAPPING MODE 



The electronic Raman scattering in superconductors 
are well described in Ref. |3C[| . Therefore, here we give a 
brief summary of the Raman spectra. The Raman spec- 



tra Si in superconductors measure effective density fluc- 
tuations, and its strength of the scattering is determined 
by the Raman vertex, ji. One can select the Raman ver- 
tex which allows for different projections on the Fermi 
surface. The intensity of the each Raman mode pro- 
vides information on the gap structure along the Fermi 
surface. The Raman spectra in superconductors is deter- 
mined from 



S?(u>,q-0)=lm 



(7»,7») 



<7i,l) 
(1,1) 



(4) 



Here we use the following notational convenience, 



(A,B)=TJ2Y1 Tr ^P3G(p, u n )B P3 G(p-q, 



n p 



(5) 

where the single particle Green's function, G(iw„,k), in 
the Nambu space is given by 



G~ 1 (iuj n , k) 



£ k p 3 - A(k ■ p)cn. 



(6) 



Here pi and Oi are Pauli matrices acting on the particle- 
hole and spin space, respectively. uj n = (2n+ l)irT is the 
fermionic Matsubara frequency, and £k — (fc^+£^)/2m — 
p, where p is the chemical potential. The Raman vertics, 
Bi g and B 2g are written as 



-fBlg = \/2cOs(2(/>) 

lB2g = V2sin(2</>), 



(7) 



where <fi is the angle of the wave- vector k on the Fermi 
surface. The second term of Eq. ^ is a back-flow term 
due to the charge conservation, which can be seen in the 
limit of ji — const. In this case, Raman intensity should 
vanishes because there is no density fluctuation in the 
homogeneous limit of q — > in a superconductor. 

The clapping mode is the fluctuation of the order pa- 
rameter which can be written as SAp 3 ~ exp (±2i0)(cr 1 ± 
i a 2)Pz in the Nambu space, As it is indicated in the fluc- 
tuation of the order parameter, <5A oc exp(±2i0), this 
mode can directly couple to B\ g and B 2g channels of the 
Raman spectra. The clapping mode makes the additional 
contribution to the Raman intensity given by 



S, = S? + Im 



( 7 i,<SA)(<SA,7«) 
g-i-{5A,5A) 



(8) 



where g is the coupling constant which mediates super- 
conducting state, and the coupling between the fluctua- 
tion of the order parameter and the light scattering with 
Raman vertex, 7j. On the other hand, the A\ g mode 
does not couple to the clapping mode. Therefore <SUi s 
is given by the first term in Eq. [§. Note that Sf is the 
same for all three modes (Bi g , B 2g , and Ai g ) for each or- 
der parameter. This indicates the existence of the axial 
symmetry and horizontal nodes in the order parameter, 
if there is any. 
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III. RAMAN SPECTRA OF B lg AND B 2g IN 
P-WAVE SUPERCONDUCTOR 

The order parameter of Eq. [j] has a full gap on the 
Fermi surface, therefore the bare Raman intensity in p- 
wave superconductor is same as that of s-wave supercon- 
ductor. 



S° Blg = Im(cos(2<£),cos(2c£)) = 



2tt7V(0)A 2 



4A : 



=t9(a, 2 - 4A 2 ). 
(9) 

Note that the Raman intensity is zero for the frequency 
u> < 2A due to the presence of the full gap. 

The coupling to the collective mode leads to the addi- 
tional contribution to the Raman spectra. While the A\ g 
mode has the same frequency dependences as B\ g and 
Big as far as the bare Raman intensity, S° is concerned, 
this mode does not couple to the clapping mode. The 
Raman intensity due to the clapping mode can be ob- 
tained by computing the correlation functions presented 
in the Appendix. Here we summarize the result of the 
correlation functions. 



Re(6A, lBlg ) = N(0)[uAf], 



Re(SA,SA) 



where / is given by 



dE 



N(0) 



tanh(£:/2T) 



T 



^/E 2 - A 2 4E 2 -uj 2 



(10) 



(11) 



Since the imaginary part of (5 A, cos(2c/>)) and (SA,SA) 
are zero for the frequency, to < 2A, the clapping mode 
appears as a resonance in the Raman spectra for B\ g , 
which is given by 



S^(ij) = Im 



(■j Bls ,SA)(6A, lBla 
g-i-(5A,8A) 



= 2nN(0)u 2 A 2 f(uj,T) 8{lj z - 2A 2 ).(12) 

The Raman intensity, B\ g (ui) is plotted in Fig. 2 with a 
finite impurity scattering, T = 0.1A. It is important to 
note that there is a resonance at the frequency of u> = 
V2A with a Lorentzian shape. The same is true for Bi g 
channel. 



IV. RAMAN SPECTRA OF B lg AND B 2g IN 
F-WAVE SUPERCONDUCTOR 



20 



15 



S B / N(0) 

10 



5 - 



0.4 



0.8 



1.2 

CO / A 



1.6 



2.4 



FIG. 2: Raman intensity for p-wave with a finite impurity 
scattering, T — 0.1 A. The solid line is the contribution from 
the clapping mode, the dotted line is the bare Raman intensity 



where K and E are complete Elliptical integral of first 
and second kinds, respectively. 

The contribution from the coupling to the clapping 
mode is given by 

S Big - N(0)uj A (w2 _ w k )2 + 4(wM2 . (I 4 ) 
where 



2 r° 

D(u>,T) = - \ 



dE If ,E S ,E\ , E 

— n n-r K(-) -E(-) tanh — 

AE 2 - oj 2 A \ K A' y A J J 2T 

(15) 

The frequency of the clapping mode and its damping are 
obtained as 



= 1.02 A 
J - 0.57A, 



(16) 



where the 

Im(<5A,(5A) 



r c is computed 
and Im(<5A,(5A) 



from the formula, 
is given in the Ap- 
pendix, Eq. A9. The Raman intensity due the clapping 
mode is shown in Fig. 3. It is important to note that 
the shape of the Raman intensity is not Lorentzian, 
but it has asymmetry. The Raman spectra for the 
another proposed f-wave, Eq |3| are identical to the result 
presented here for the order parameter of Eq. 0. 



We consider the f-wave order parameter of Eq. |2|, and 
the cylindrical Fermi surface independent of k z . Since 
there is a node along the k z direction, we expect the 
finite bare Raman intensity, (7^,7^) at low frequency, 



r.O 



- Im(cos {2(f>), cos (2(f))) 
AN(0)A 



2 A' 



tanh 



4T' 



(13) 



CONCLUSION 



Within the framework of the weak-coupling BCS 
theory, we have studied the Raman spectra of two- 
dimensional p-wave (Eq. 0) and f- wave (Eqs. and |^) 
superconductors; the candidate for S^RuCi. The clap- 
ping mode with angular momentum ±2, parallel to the 
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APPENDIX A 



1. the clapping mode and its contribution on 
Raman, £>i 9 ,2 g in p-wave superconductor 
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FIG. 3: Raman intensity for f-wave. The solid line is the 
contribution from the clapping mode, the dotted line is the 
bare Raman intensity, and the dashed line is the total Raman 
intensity 



c-axis is common to these ground states, which break 
the chirial symmetry. We have shown that the clapping 
mode has the frequency of ^/2A(T), and 1.02A(T) for p- 
wave and f-wave superconductors, respectively. We have 
also shown that the clapping mode couples both sound 
waveQ, and B\ g and Bi g modes of the Raman spectra. 
While the coupling to the sound wave is very small, the 
clapping mode appears as a huge resonance in Raman 
spectra. Therefore, if the present experimental difficulty 
is overcome, the Raman spectroscopy provides a unique 
window to probe the clapping mode. 

Investigating the clapping mode in Raman spectra will 
also probe the order parameter of Sr2Ru04. The sin- 
gle model of the f-wave shows only contribution from 
the clapping mode of f-wave order parameter( Fig. 3 
). On the other hand, the multigap model proposed 
by Zhitomirsky and Rice should be the weighted sum 
of the p-wave (Fig. 2) and f-wave (Fig. 3) contribu- 
tion. It is important to note that the maximum gap of 
A zr (0.2 ~ 0.5)Abcs on the a- and /3-bands, while 
the maximum of A on 7-band with p-wave is set to be 
Ascs which is the energy scale we used in our figures. 
Therefore, for the multigap model, the contribution from 
f-wave order parameter, Eq. |^ should be peaked around 
(0.2 ~ 0.5)A (with the contribution from p-wave around 
y/2A), while single gap model gives the peak around 
1.02A. The Raman spectroscopy can discriminate the 
multigap model from the single gap model. 
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The fluctuation of the order parameter corresponds to 
the clapping mode, and it can be written as SAp 3 ~ 
e ±2t ^(ai ± i(J2)p3- Its coupling to the 73 i g mode of Ra- 
man spectra, and its propagator, (5 A, 5 A) are obtained 
by computing the following correlation function within 
weak coupling theory. 

(*A,cos(20))(iw„,q) = T^^Tr[ < 5A /O3 G(p,w„)cos(20) 
(5A,SA)(iu v ,q) 



xp 3 G(p - q, iiv n 
T^2J2Tr[5Ap s G(p,cj n )SA 



xp 3 G(p - q, iuj n 



(Al) 



After summing over the frequency, the correlation 
functions are written as 



Re(5A, 7mg ) = N(0) 



dE 



loA 



1 



= [N(0)wAf], 



E 

V £ 2 -A2 4E*-u* 2T 
(A2) 



where f(ui,T) is given by Eq. [n|. There is no imag- 
inary part of the above correlation function. On the 
other hand, the correlation function of the order param- 
eter fluctuation is given by 



Re(5A,8A) = N(Q) 



OO 



dE 2E 2 



A 

UJ 2 

— L 

2 



A 2 , E 

— s ?r tanh — 

A 2 UJ 2 - 4£ 2 2T 



f 



where 



N(0) 



A 



\JE 



dE , E 

— tanh —= 

A 2 2T 



(A3) 



(A4) 



Now the imaginary part of the above correlation function 
is 



r 00 

Im(SA,SA) = -N(0) — / dE 
2w J A 



IE 2 -A 2 , E 
tanh — 



VE 2 - A 2 -T 
x (S(ui - 2E) + 5{u + 2E)) 



N(0) 



7T LO 2 - 2A' 



_tanh— 0(w-2A). (A5) 



Note that Im(dA,SA) is for to < 2A, which is the 
regime we are interested in. Therefore, the Raman in- 
tensity for Big as well as B 2g for p-wave for u < 2A 
should be 



Im 



Im- 



2N(0)uj 2 A 2 f 



/ ([ 7 ,<5A P3 ])([(5Ap3,7]) 

\g-i - {{5Ap 3 ,SAp 3 }) J ^-{lo 2 -2A 2 )+iT 
= 2TTN(0)uj 2 A 2 fS(uj 2 - 2A 2 ),(A6) 
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where a finite T could come from the impurity scattering. 
The clapping mode appears as a huge resonance with 
Lorentzian shape (with a finite T) at the frequency of 
\/2A shown in Fig. 2. 



2. the clapping mode and its contribution on 
Raman, Bi 9 ,2 9 in f-wave superconductor 



Using the correlation functions, Eq. Al, we obtained 



the following results for the f-wave superconductor. 

2iV(0)wA f°° dE 



Re(5 A, cos (20)) 



o u 



4E 2 



x (A^tanh^ (A7) 



Using the definition of D(lo,T), Eq. |T^, the above corre- 
lation function can be written as 



Re(8A, cos (20)} = N(0) [uAD(w, T)] 



(A8) 



On the other hand, the propagator of the clapping mode 
is given by 



Re(SA,6A) = g- 1 - 



2N(0) 



dE 



AE 2 

,E 



[ —A(^)- A 2 B(^-) )tanh — 
2 V A ; V A 7 / 2T 



Im(SA, SA) 



g- 1 - N(0) 
N(0) 



LO 2 - CJ 2 



D(u,T) 



LU 2 A (^-)-2A 2 B(^ 
K 2A' v 2A 



LO 

x tanh — — , 
4T' 



(A9) 



where 



A(k) ee -( K (k)-E(k)) 
B(k) 



A 
1 

3A 

2A 2 B(k)/A{k). 



[{k 2 + 2)K{k) - 2(k 2 + l)E(k)] 



(A10) 



Therefore, one can determine the position of the clapping 
mode, ujf}, using the following relations. 



2k 2 = - 2 



k 2 K(k) - 2k 2 E(k)\ 
K{k) - E{k) ) ' 



(AH) 



where k = u>^/(2A). This gives, 



? = 1.02A. 



(A12) 
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